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on : We propose a formula for flux of extremely high energy cosmic rays (EHECR) through 

decay of superheavy particles. It is shown that EHECR spectrum reported by AGASA is 
reproduced by the formula. The presence of EHECR suggests, according to this approach, 
the existence of superheavy particles with mass of about 7 x 10 11 GeV and the lifetime of 
about 10 9 years. Possibility to obtain a knowledge of f2 of the universe from the spectrum 
of EHECR is also pointed out. 
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• The unexpected energy spectrum of the extremely high energy cosmic rays (EHECR) with 

the energy above 10 9 eV has been reported by Akeno Giant Air Shower Array (AGASA) 
collaboration which used the updated data set[Q]. The existence of EHECR has been known 
for about 30 years ||, among which the highest energy of EHECR, (3.2 ± 0.9) x 10 20 eV, was 
recorded in Fly's Eye ||. When we regard nucleons or nuclei as constituents of EHECR, the 
attenuation length is estimated less than ~ 100 Mpc due to Greisen-Zatsepin-Kuz'min (GZK) 
effect j|, [|. Thus if distances from Earth to any sources of EHECR are over 100 Mpc, it is 
difficult to explain EHECR in terms of nucleons or nuclei. One may expect that sources of 
EHECR would exist within about 100 Mpc from Earth, but such sources have not been found. 
Here we can make a brief list of yet unsolved problems on EHECR as follows: (1) the chemical 
compositions of EHECR are unknown, (2) the sources of EHECR are unclarified, and (3) the 
shape of energy spectrum reported by AGASA is unexplained. 

In order to solve these problems, several ideas have been proposed. There exist mainly two 
approaches based on astrophysical aspect and particle physical aspect. From the astrophysical 
^ ■ aspect, number of acceleration mechanisms of the ordinary particles have been proposed. For 

example, the protons are accelerated to the extremely high energy by the relativistic jets from 
AGN H . The production of EHECR by the Gamma- Ray Bursts is also considered Q • There 
are, however, no particular astronomical objects in the directions of EHECR [|], [|. 

On the other hand, from the particle physical aspect, the production mechanism of EHECR 
has been proposed |§. A basic idea is that EHECR are produced by decay of superheavy 
particles. In this case, EHECR are consisted of either the standard particles (proton, photon 
or neutrino) or the new particles (for example, neutralino or gluino). Furthermore, instead of 
decay products of superheavy particles, quasi-stable superheavy particles can also be considered 
as constituent of EHECR in itself. For example, colored monopoles have been examined in 
ref.§. 

When decaying superheavy particles are considered as sources of EHECR, it is plausible to 
imagine the mass of the superheavy particles larger than about 10 12 GeV, since the highest energy 
of EHECR is (3.2± 0.9) x 10 20 eV. There are many models in particle physics which introduce such 
superheavy particles, for example, supersymmetric grand unified theory or the see-saw model 
for neutrino mass. However, the superheavy particles that should decay into EHECR have to 
have a long lifetime which is nearly or greater than the present age of the universe, because, if 
the superheavy particles decay fast, the remnants of the superheavy particles are very little, and 
as a result, we cannot observe EHECR produced via the decay. Although the order of lifetime 
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is usually proportional to the inverse of mass of the particle, some kind of models |Uj] permits 
quasi-stable superheavy particles contrary to their heaviness. 

It is attractive to make clear of the relation between the EHECR problem and the particle 
physical aspect, because we can expect that the feature of EHECR leads us to ones beyond the 
standard model of particle physics with a large mass scale larger than about 10 12 GeV. 

In this paper, we examine the problem of EHECR from the particle physical aspect, although 
we do not reject at the present time a possibility to explain EHECR from astronomical origins 
like quasisteller objects which are recently considered seriously in ref. [pj|| . We propose a proper 
formula for the flux of EHECR in the beginning. This formula should be used when EHECR 
are produced via the decay of the superheavy particles. Reproducing the spectrum of EHECR 
observed at AGASA by using our formula, we show that the lifetime of the superheavy particles 
is suggested to be about 10% of the Hubble time when their mass is 7 x 10 11 GeV. Furthermore, 
it is also pointed out that the possibility to obtain a knowledge of the omega parameter, £1$, of 
the universe from the energy spectrum of EHECR. 

At first, we show the formula for the flux to the cosmic rays as decay products of superheavy 
particle X. Let / be the particle in decay products of X and produce Extensive Air Showers 
(EAS) when it reaches to the atmosphere. We denote the number density of the source at time 
t e as nx(t e ) and the partial decay constant for a decay mode including / as Tf. Then the 
production rate of / at time t e is determined by Tfnx(t e )- The general expression of the flux 
of cosmic rays, after taking average about the angular distribution, is given in the Robertson- 
Walker metric with the scale factor a(t) as follows: 

, . 1 f*o , fa(t e )\ 3 , ,dT f dE e 

J(E obs ) = — dt e nx (t e )—L—Z-, 1 

4?r Jtmin \a{t )J dE e dE obs 

where E obs is the observed primary energy of cosmic rays, and dT j /dE e is the energy distribution 
of cosmic rays emitted by the source at time t e . The present time is represented by to, while 
t m i n , is determined by a physical condition that to — tmin is smaller than an attenuation time of 
/• 

The number density of the superheavy particles nx(t) is evolved by the following Boltzman 
equation, 

*^ = -r wl .)' w («], (2) 

where Tx is total decay width of X. This equation gives 

a(t) 

By substituting eq. @ into eq. @), we obtain 

J(E obs ) = dt e ex V [r x (t -t e )]-f—^- . (4) 

47r Jt min 0-E e dti obs 

Here we would like to make a comment on this formula. In the previous works, similar formulae 



nx(t) = (^\ nx(t )exp[r x (t -t)]. (3) 



have been used ||, 12]; however, the effect of the decrease in the number density of parent 
particles has been neglected, since their lifetime has been assumed to be greater than the age 
of the universe. On the other hand, in eq. (||), this effect is included in the form of exponential 
damping factor, exp[Tx(to — t)], since X is supposed to have the same order lifetime as the 
present age of the universe. We should use eq. (Q) rather than the formula used in previous 
works when we estimate the flux of EHECR produced by the decay of superheavy particles. 
This is because the number density of parent particles decreases due to their decay. This effect 
is very important to reproduce the spectrum of EHECR as shown later in this paper. 
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Our formula, eq. (||), can be applied to the case that / can travel almost freely through the 
cosmic background radiation, though the attenuation due to the GZK effect has been neglected. 
For example, we can take / to be neutrino or neutralino. The difference between them is the 
translation rate R C 2a of EHECR to EAS when they reach at Earth and produce EAS. 



R, 



number of EHE air shower events 



c2a 



number of EHECR incident on Earth 
Then the flux of EAS, J(E obs )\EAS, is derived from the flux of EHECR, J(E bs)\EHECR, as 



J(E, 



obs)\EAS 



J(E, 



obs ) I EHECR Rc2a 



(5) 



(6) 



Neutrino gives the translation rate larger value as an order of R C 2 a = 10 -6 |l3| than neutralino. 

Now we consider a case that EHECR are taken to be neutrinos which are produced by the 
two-body decay of X. Its mass, Mx, can be expected about twice the highest energy of EHECR. 
Since the highest energy recorded by Fly's Eye is (3.2± 0.9) x 10 20 eV, we take Mx as 7x 10 11 GeV 
in the following. The energy spectrum right at the emitted point is to be the monochromatic 
one in the two-body decay, 



dr 



/ 



dE P 



T f 5{E e -M x /2). 



(7) 



Hereafter, for simplicity, we limit ourselves to consider such a case that there is unique 
channel for two-body decay, i.e. X — > u+ some particle, and then Tj becomes just the total 
decay width T tot . Since the kinematics tells that the energy spectrum is monochromatic in the 
present case, we can obtain the energy distribution of EHECR without the detailed description 
of the interaction. Thus we will be able to discuss the problem of EHECR rather in the model- 
independent way. 

The red-shift relation between E e and E a b s is represented as (1 + z)E Q bs = E e . Then the 
monochromatic energy condition becomes 



6{E e - M x /2) 
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E obs \df(t)/dt\ 



5(t-t a )6(z), 



(8) 



where the step function 9{z) is necessary to satisfy the condition z > 0. In eq. ((^), f(t) = l+z(t), 
and t a is defined as a solution of f(t a ) = Mx /2E t, s . The expression of f(t) is different for each 
of three cases, open universe of f^o < 1, flat universe of f^o = 1 and closed universe of £Iq > 1. 
By using eqs. (Q) ~ (^), we can calculate the energy spectrum of EHECR. For the moment, we 
consider the case for Qq = 1, where f(t) = (to/t) 2 / 3 . In this case, the spectrum of EAS is given 
by 



J (E bs) 



. V2T X VEoVs \ 2Y X 



2E, 



obs 



M x 



3/2 



(9) 



The AGASA's data above the GZK cutoff jjj distributes over the energy region from 10 19 - 8 eV to 
10 20 ' 5 eV, which corresponds to the change of 1 to 5 for 1 + z. This means that the most remote 
position of X is about 7 Gpc in flat universe case. This distance is certainly shorter than the 
mean free length of supposed neutrinos on the collisions with background photons or neutrinos 
[14]. Thus we can practically put t m i n into in eq. (||). We see from eq. (||) that the spectrum 
of EHECR is determined by the three parameters Mx, r = Tx/Hq, and R C 2a n x(to)- Since R C 2a 
appears as a product with nx(to) in eq. @, fortunately an ambiguity of R C 2a does not affect 
our estimate of the spectrum. When we reproduce the spectrum reported by AGASA, we must 
multiply the factor E^ bs to eq. @. 
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Fig. |l] displays the curves of log 10 ( J(E f, s )E^ bs ) for the case that EHECR are to be neutrinos 
produced by the two-body decay of X. The flux J(E b s ) is calculated by eq. (^). The end point 
of the energy spectrum locates at Mj/2. The energy of EHECR produced at time t a is lowered 
to E e /(l + z) by the red shift. The peak around 10 20 eV is reproduced by synergistic effect 
between the decreasing factor of X, exp[Tx(to — and the expanding factor of the universe, 
(io/ii) 2 / 3 in eq. (|j). We see remarkable conformity of the energy spectrum of EHECR by eq. (^) 
with the one reported by AGASA. 

The location of peak of the spectrum shown in Fig. [I] is derived from eq. (||) as 

^(peak) = ^^) 2/3 . (10) 

In order that the peak appears in Fig. ffl, r > 3.5 is needed. We see from the data point given 
by AGASA that the peak seems to locate at E b s ~ 10 20 eV, though this location is vague in the 
data. When the inequalities 20.1 < log 10 E Q b s < 20.4 are imposed, we obtain 

3.5 (^HT 11 - 4 ) 1 ' 5 < r < 3.5 (^lO- 111 ) 1 ' 5 (11) 

where Mx is given in the unit of GeV. For Mx = 7 x 10 11 GeV, we find 5.8 < r < 16. The 
inequalities turn the lifetime of X into tx = (0.09 ~ 0.3) to- 

When we see AGASA's data [|lj in detail, there is no event in the energy bin between 10 20 ' 2 
and 10 20 3 eV. We anticipate, however, appearing of one event at this energy bin from our present 
consideration. But if there remains this vacancy in future, this spectrum should be interpreted 
from our standpoint as a suggestion that two kinds of superheavy particles with different masses 
exist. The spectrum below 10 20,2 eV is due to superheavy particles with mass of 3 x 10 11 GeV, 
while the one over 10 20,3 eV is due to another superheavy particles with mass of 7 x 10 11 GeV. 

In Fig. ||, the curves of log 10 ( J{E ^ s )E^ )bs ) are shown for J7 = 0-5, 1, and 2. We put r = 10 
for Mx = 7 x 10 11 GeV as an example. Statistics of the current data above 10 19,8 eV is not 
enough to derive information of Oo- We hope from this point that the statistics of the EHECR 
observation increases in future. 

So far we have considered two-body decay case. Here we address a question how the en- 
ergy spectrum changes for multi-body decay case. As an example, we show the curves of 
log 10 ( J(E f )S )E^ bs ) for three-body decay due to four-Fermi interaction. We can see from Fig. || 
that the curves for three-body decay case do not reproduce the data so well as those for two-body 
decay case in Fig. [j]. 

In summary, we have examined the problem of EHECR from the particle physical aspect, 
proposing a suitable formula to the flux of the cosmic rays through decay of superheavy particles. 
Then the energy spectrum of EHECR has been reproduced satisfactorily as that of the neutrinos 
via two-body decay of X with mass of 7 x 10 11 GeV and lifetime of around 0.1 R C 2a^x(to) ~ 
10~ 14 has been required to reproduce the energy spectrum of EHECR, where Q,x = Px/ Per as 
is usually defined. Our analysis is not confined exclusively to the neutrino case as mentioned 
before. We can take neutralino instead of neutrino. The difference between neutralino and 
neutrino is only the value of R C 2 a , and Rica becomes smaller for neutralino. Thus, &x(to) turns 
out to be much larger for neutralino case than 10 -8 for neutrino case. Furthermore, we have the 
possibility to derive a knowledge of the omega parameter, f^o, of the universe from the energy 
spectrum of EHECR. We expect to acquire further information of superheavy particles as well 
as Slo °f the universe from future observations where more data will be accumulated to give 
higher statistics. From this aspect, future detectors like the HiRes, the Telescope Array, the 



Pierre Auger, and the Owl [15] are greatly awaited. 

One of the authors (Y. C~) thanks Tetsuya Hara for his kind correspondence. 
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Figure 1: Curves of J(E f, s )E^ bs for the two-body decay case. Data points and error bars are 
taken from AGASA data given in Q. Input parameters are taken as follows; Mx = 7x 10 11 GeV, 
r = r x /H = 9.5, 10, and 10.5, ft = 1, Px{t ) = M x n x (t ) = 10"Vr(to), and R c2a = lO" 6 . 



6 




19.5 20.0 20.5 

Log 10 (E obs eV) 



21.0 



Figure 2: Curves of J(E b s )E^ bs for the two-body decay case. We take the omega parameters 
as Qq = 0.5, 1, 2, and r = 10. Other input parameters are taken to be the same as in Fig. [l} 
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Figure 3: Curves of J(E i, s )E^ bs for the three-body decay case. We take the ratio of Tx/Hq as 
r = 10, 10.5, and 11.5. Other input parameters are taken to be the same as in Fig. |]. 
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